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The current induced local spin polarization due to weak Rashba spin-orbit coupling in narrow 
strip is studied. In the presence of longitudinal charge current, local spin polarizations appear in 
the sample. The spin polarization perpendicular to the plane has opposite sign near the two edges. 
The in-plane spin polarization in the direction perpendicular to the sample edges also appears, but 
does not change sign across the sample. From our scaling analysis based on increasing the strip 
width, the out-of-plane spin polarization is important mainly in a system of mesoscopic size, and 
thus appears not to be associated with the spin-Hall effect in bulk samples. 



In a spin-orbit coupled electron system, an external 
electric field can induce a transverse spin current, giv- 
ing rise to the so-called spin Hall effect (SHE). The SHE 
may offer a new way to control electron spins in semicon- 
ductors, and so have potential applications in spintronic 
devices. Depending on its origin, the SHE is generally 
divided into two categories: the extrinsic SHE, which 
originates from spin-dependent electron anomalous scat- 
tering by impurities, and the intrinsic SHE, which oc- 
curs even in the absence of impurities. The extrinsic 
SHE, was first proposed by D'yakonov and V. I. Perel 1 
in 1971 and reexamined recently by Hirsch 2 and Zhang 3 . 
The intrinsic SHE was predicted by Murakami, Nagaosa, 
and Zhang 4 for p-type semiconductors and by Sinova et 
al£ for n-type semiconductors in two-dimensional het- 
erostructures. The intrinsic SHE has attracted much 
theoretical interesta 6 ! 7 ! 8 ! 9 ! 10 ! 11 . 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 . Very re- 
cently, two independent groups have reported experimen- 
tal evidenc e) 19 ! 20 that an electric field can cause out-of- 
plane spin accumulations of opposite sign on opposite 
edges of semiconductor films, which is considered to be 
a signature of the SHE. Several analytical and numer- 
ical works have been published on the subject of spin 
accumulation in a semiconductor with spin-orbit cou- 
pling. Governale and Zulicke2i were the first to investi- 
gate spin accumulation. They studied the spin structure 
of electron states in a quantum wire with parabolic con- 
fining potential and strong Rashba spin-orbit coupling. 
Usaj and Balseiro 2 ^ showed that in a semi-infinite sys- 
tem with spin-orbit coupling, a current flowing parallel 
to the edge induces a net magnetization close to the edge. 
Using the Landauer-Biittik formula for a tight-binding 
model, Nikolic et al~ showed numerically that in a two- 
dimensional bar with a width of 30 lattice constant, the 
Rashba spin-orbit coupling can induce opposite spin ac- 
cumulation near the two edges, which is qualitatively sim- 
ilar to that observed in the experiment. In order to clarify 
whether such spin polarizations are related to the SHE, 
it is important to investigate their scaling behavior with 
increasing sample size, and to reveal the parameters that 
control the relative magnitude of the spin accumulation 
or polarization. 

In this paper, the electron wave function in a continu- 



ous model is obtained for an infinite long conducting strip 
with finite width L. Using the Kubo formula, we show 
that a longitudinal electrical current induces both out-of- 
plane spin polarization (S z ) and in-plane spin polariza- 
tion (S y ). Near the two edges, the spin polarization S z 
has opposite sign, whereas S y has the same sign. When 
sample width L increases, its scaling behavior indicates 
that S z near the edges decreases and S y becomes dom- 
inant for given fixed electrical current density. There- 
fore, the out-of-plane spin polarization is an effect due 
to boundary reflections from the two opposite edges, and 
appears not to be related to the SHE in a bulk sample. 

Let us consider a system of a two-dimensional (2-D) 
infinite long conducting strip with finite width L. The 
Hamiltonian for the system with Rashba spin-orbit cou- 
pling can be written as by 

k 2 

H = — + X(a x k y - (T y k x ), (1) 

where A is the coupling constant of spin-orbit interac- 
tion, a x and a y are the Pauli matrices, m is the electron 
effective mass, and wc take units with Ti — 1. 
The eigenstates of plane waves are 

\E T ,k>=\E T ,k x ,k y >=— e ife ^ 1 j, (2) 

where <f> — arctan(fcj,/fc x ), + (— ) labels lower (higher) 
energy eigenstate with eigenvalue E± — k 2 /2m =F Afc for 
a given k. 

Assuming hard-wall boundary conditions, the wave 
function at the two edges (y = and y = L) is zero. 
Since the system is uniform along the x direction, k x 
commutes with the Hamiltonian and is a good quantum 
number. We can write a eigenstate, with eigen-energy 
E, of the system as a superposition of four plane waves, 
with same E and k x . Suppose the system is in universal 
region as defined in Ref. [5], the wave function near the 
Fermi level is given by 

*(E,k s ,y) = \E,k x > 

= a k JE,k x ,ky > +p k JE,k x ,-k- > (3) 
+ j k JE,k x ,k+ > +5 kx \E,k x ,-k+ >, 
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where k^ — \J k ±2 — k 2 and = ±Ato+\/A 2 to 2 + 2mE 
with boundary conditions ^(k x ,0) = ^(k x ,L) = 0. One 
can solve the boundary conditions and find the eigen- 
values of k x , which are a discrete set of values in the 
interval of (—kp, kp), here kp = k + with E = Ep (the 
Fermi energy). In Usaj and Balseiro's work, there is only 
one edge, the eigenfunctions are propagating waves writ- 
ten as a superposition of one incident and two reflected 
waves. k x can take any value between (—kp,kp). In 
our current study, the interference due to the two edges 
of the strip limits number of eigenvalues for k x at Fermi 
level, which could inject rather different physics for the 
problem. 

While the four plane waves have different spin polariza- 
tions within the two-dimensional plane, their interference 
leads to nonzero out-of-plane local spin density. The lo- 
cal spin polarization depends on the sign of the conserved 
longitudinal wave vector k x . For any given energy E and 
a positive eigenvalue kEnx for k x , — kp n x is also an eigen- 
value for k x . In the ground state where both positive and 
negative k x states are occupied, the total local spin den- 
sity is zero since the contribution of each spin band is 
zero. However, if there are a longitudinal current flow- 



ing in the strip, which causes a small shift of the Fermi 
circles. The numbers of occupied states with positive k x 
and negative k x are no longer equal, which can induce 
net spin polarizations in the strip. 

The net local spin polarization can be calculated using 
Kubo formula 2 ^, 
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{< E',k x \\a{y)\E,k x >< E,k x \v x \E',k x >} 
(E — E') (E — E' — iS) 



(4) 



where £ is the electric field, and v x is given by v x = 
k x /m — \a y . Here, we wish to point out that the con- 
tribution to the spin polarization comes from the Fermi 
surface, in contrast to the intrinsic SHE, which comes 
from the contribution of all occupied states. It is easy to 
find that the spin polarization given by Eq I0J diverges in 
the clean limit. To overcome this problem, we consider 
a constant longtidudinal current density I/L is driven 
through strip. We calculate the ratio between the spin 
polarization and the current density 
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which is a finite quantity. 

All the coefficients in Eq.[3|can be determined numeri- 
cally. However, we found that the eigenfunctions (stand- 
ing waves) can not be expressed as a superposition of 
the two eigenfunctions^ of same k x obtained in the case 
of only one edge. We also found that, for \k x \ < kp, 
\a\ = \0\ and | 7 | = \S\; for k F + > \k x \ > kp, | 7 | = |<S|. 

Plots in Fig. show local spin polarizations s z = | < 
Ef, k x \a z \E F , k x > and s y = | < E F ,k x \a y \E F ,k x > 
as functions of position y in all the eigenstates that have 
positive k x with eigen-energy at the Fermi level for three 
different values of Rashba coupling A, which is in the 
units of k F . The width of the strip is set to be 8/kp, 
where kp is the Fermi wave vector when there is no 
spin-orbit coupling, kp is related to electron density in 
the sample by kp = 2im. Using typical value of two- 
dimensional electron density 10 12 cm -2 [see Ref. 26], we 
estimate kp 10 8 /m and L is around 80 nm. The spin 
polarizations s z and s y vanish at the two edges as re- 
quired by the boundary conditions. For each eigenstate, 
Sz{y) = s z (L - y) , whereas s v (y) = s y (L - y). We 
have also obtained s x and it is zero across the sample. 
Without spin-orbit coupling, k y is quantized to values 

For each k yn , the 



mr/L, where n = 1,2,3, 



eigenstates for two spin directions are degenerate. When 
spin-orbit coupling is in presence, the two spin bands 



are no longer degenerate. However, the spin polariza- 
tion increases as the Rashba coupling increases. When 
we further increase A, some values of k x are larger than 
Fermi wave vector of the higher spin band, as shown in 
Figs.^e) and (f). Under this case, decaying waves show 
up in the wave functions along the y-direction for the 
higher spin band. 

The net spin polarizations are calculated by using the 
Kubo formula in Eq. (0J. Figures |2Ia), (c) and (e) show 
the net S z (y)L/I when we sum the contribution from all 
the positive k x modes at the Fermi level. The longitudi- 
nal charge current also induces a local in-plane polariza- 
tion Sy&Z^, as shown in Figs. Etb), (d) and (f), whereas 
S x = 0. Unlike S z , S y has the same sign across the sam- 
ple. At weak Rashba couplings A = 0.01 or 0.05, we see 
from Figs. |2 (a)-(d) that S z is one or two order greater 
than S y in magnitude. With increasing the Rashba cou- 
pling, S y increases much faster than S z . As a conse- 
quence, S y becomes comparable to S z at relatively large 
Rashba coupling A = 0.1, as seen from Figs.|2Ie) and (f). 

We have also obtained the local spin polarization for 
the case of L = 16 /kp. Figure |21 shows the results for 
A =0.01, 0.05, and 0.1. There are 10 eigenvalues for k x 
for each A. Similarly to Fig. [21 the magnitude of spin 
polarization S Z L/I near the edges increase as A increases 
when A is small. However, when A is large (A = 0.1), 
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FIG. 1: s z [figure (a), (c), and (e)] and s y [figure (b), (d), (f)] 
as function of position y for eigenvalues of positive k x at the 
Fermi level. s z and s y are in units of |. L — 8/kp, where &f 
is the Fermi wave vector when there is no spin-orbit coupling. 
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FIG. 2: Total local S Z L/I and S y L/I as a function of y for 
L = 8/fci?. In this and following figures they are in units of 
h/(2ek%) In (a) and (b) A = 0.01. In (c) and (d), A = 0.05. 
In (e) and (f), A = 0.1. 



magnitude of the polarization S z near the edges becomes 
smaller and large oscillations appears deep inside of the 
sample. S y also increases as A increases. When A = 0.1, 
in-plane spin polarization S y dominates out-of-plane po- 
larization S z . Comparing with the results of L — 8/kp, 
S Z L/I is larger near the edges for A =0.01 and 0.05. But 
when A =0.1, it is smaller. S y L/I is larger for all A's 
we chose. We conclude that for fixed finite sample size, 
the out-of-plane spin polarization S z dominates at rela- 
tively weak Rashba coupling and in-plane spin polariza- 
tion overwhelms for relatively strong Rashba coupling. 
It is interesting to examine how the spin polarization 
changes as the sample width L increases. In Fig. 0Ja), 
we show S Z L/I for various L with A = 0.05. The peak 
magnitude of S Z L/I near the y=0 edge increases as L 
increases at small L and decreases as L increases at large 
L. The plots for A = 0.1 in Fig. 0Jb) show a similar pat- 
tern, but the width, which has the biggest S Z L/I near 
the edge, is shorter than that of A = 0.05. 

In summary, we showed that in a 2-D narrow semicon- 
ductor strip with weak Rashba spin-orbit coupling, local 
spin polarization could be induced by a steady longitudi- 
nal current, and it is originated from the wave functions 
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of the electrons at the Fermi level. The charge current 
along the strip induces both out of plane and in plane 
local spin polarizations. Near the two edges, the spin po- 
larization S z has opposite sign, whereas S y has the same 
sign. When the sample width L increases, the peak mag- 
nitude of S Z L/I near the edges increases at small L and 
decreases at large L for weak A. And at large L, our nu- 
merical results indicate that S y L/I becomes dominant. 
From our scaling analysis based on varying L, the out-of- 
plane spin polarization is important mainly in systems of 
mesoscopic sizes, and thus appears not to be associated 
with the SHE in bulk samples. 
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FIG. 3: Total local S Z L/I and S y L/I as functions of y for 
L = 16/fc F . In (a), A = 0.01. In (b), A = 0.05. In (c), 
A = 0.1. 
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FIG. 4: Plots of S z {y)L/I for various L with A = 0.05 [panel 
(a)] and 0.1 [panel (b)]. The value of L of each curve can 
be identified from the right end of the curve. In panel (a), 
£=4,8,12,16,24; in panel (b), £=4,8,12,16,20. 



